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I use an interpolation formula, introduced recently by Guerra and Toninelli in order to prove the 
existence of the free energy of the Sherrington-Kirkpatrick spin glass model in the infinite volume 
limit, to investigate numerically the finite size corrections to the free energy of this model. The 
results are compatible with a (1/127V) \n(N/No) behavior at T c , as predicted by Parisi, Ritort and 
Slanina, and a behavior below T c . 

PACS numbers: PACS numbers: 75.50.Lk, 75.10.Nr, 75.40.Gb 



Many years after their experimental discovery, spin 
glasses remain a challenge for experimentalists, theoreti- 
cians and more recently computer scientists and mathe- 
maticians. Numerical simulations have been used heavily 
in order to investigate their physical properties. Numer- 
ical simulations are obviously limited to finite systems. 
Simulations of spin glasses are indeed limited to very 
small systems, due to the need to repeat the simulation 
for many disorder samples (this is related, at least for 
mean field models, to the lack of self-averaging), and to 
the bad behavior, as the system size grows, of all known 
algorithms. A detailed understanding of finite size effects 
of spin glass models is accordingly highly desirable. The 
problem is also interesting in its own sakepj 0, 0. 

Here I study the finite size behavior of the Sherrington- 
Kirkpatrick modelpj (SK model), a well know infinite 
connectivity model, introduced originally in order to have 
a solvable starting point for the study of "real" finite con- 
nectivity spin glasses, and that turned out to have a com- 
plex fascinating structure, to the point of becoming^! "a 
challenge for mathematicians" . 

The partition function of the N sites Sherrington- 
Kirkpatrick model is 
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where T is the temperature, the Cj's are Ising spins, 
and the Jjj's independent, identically distributed, Gaus- 
sian random numbers with zero mean and unit square 
deviation. In the paramagnetic phase, the finite size be- 
havior of the disorder-averaged free energy f N (T) can be 
computed, using the replica method, as an expansion in 
powers of l/N, as shown by Parisi, Ritort and Slanina0. 
One starts from the equation ,6] 
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where r = (T 2 — l)/2, the field q is a real symmetric 
n x n matrix, with q a ^ a = 0. The matrix q has been 
rescaled by a factor 1/T 2 (namely q = q/T 2 ), and the 
terms of order q 5 and higher have been omitted from the 
effective Hamiltonian TL{q). In the paramagnetic phase, 
on can expand the integrand around the saddle point 
q_a, b = 0. Keeping the quadratic term only in TC, one 
obtains 
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Treating perturbatively the interaction terms in Ti one 
builds J| a loop expansion for the finite size corrections 
to the free energy. The k loops term goes like 1/N k , 
with the most diverging contribution as T — > 1 (namely 
cx \/ [N k T^ k ^)) coming from the order q 3 term in the 
Hamiltonian. Summing up these contributions, one ob- 
tains at the critical temperature 
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The computation [jj] of the constant /(-1) requires a 
non-perturbative extrapolation. The various prescrip- 
tions tried for this extrapolation gave unfortunately quite 
different values for in the approximate range 

[-0.2, +0.2]. 

It is not known how to extend the above analysis to the 
spin-glass phase below T c . Numerical works indicate that 
the ground state energy (or zero temperatures internal 
energy) scales like [7I laTfa. fiol [Tlj| ejv - cx A~ 2 / 3 (this 
result is exact for the spherical SK model 12]), like the 
internal energy at T c 0. 
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In this brief report, I introduce a numerical method to 
compute the finite size corrections to the free energy of 
the Sherrington-Kirkpatrick model, based on Guerra and 
Toninclli interpolation method. Guerra and Toninelli in- 
troduced the partition function 
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that involves a parameter t that interpolates between 
the SK model with N sites (t = 1) and a system of two 
uncoupled SK models with Ni and N 2 — N — N 1 sites 
(t = 0). In what follows 2Vi = N 2 = N/2. The Js, J's 
and J"s are independent identically distributed Gaussian 
random numbers. It is easy to show that 
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The right hand side of equation[S]can be evaluated with 
Monte Carlo simulation. I use the Parallel Tempering al- 
gorithm, with T £ [0.4, 1.3] and uniform AT = 0.025. A 
total of 2 10 5 sweeps of the algorithm was used for every 
disorder sample. The quenched couplings have a binary 
distribution in order to speed up the computer program 
(as shown in reference |l|, the leading finite size correc- 
tion is the same for the binary and Gaussian couplings). 
Systems of sizes N from 128 to 1024 have been simu- 
lated with 128 disorder samples for each system size (but 
for N — 1024, where I used 196 samples). The integra- 
tion over t was done with the trapezoidal rule, with 39 
non uniformly spaced points. Integrating with only half 
of the points makes a very small effect on the integrand 
(smaller than the estimated statistical error). 

Figure ^ shows the integrand D(t) as a function of t 
for the largest system and several temperatures. The 
integrand is concentrated around t — 0, and I have cho- 
sen the discretization of t accordingly. One notices that 
T>(t = 0) is more and more negative as T decreases, as 
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FIG. 1: (Color online) Minus T>(t) as a function of the inter- 
polation parameter t (both in logarithmic scale) for N = 1024 
and temperatures 0.4, 0.6, . . . , 1.2. 



predicted by the formula V(t = 0) = — < ((712) 2 >, and 
that V{t — 1) is weakly dependent on T, as expected 
from the identity V(t = 1) = 1/(N - 1)(< (q 12 ) 2 > - 1), 
which is weakly dependent on T (for not too small T's) 
since < ((712) 2 > is small compared to one. 

In the low T phase, a remarkable scaling is observed if 
one plots the ratio V(t)/D(t = 0) as a function of tN 2 / 3 , 
as shown in figures It means that, to a good approx- 
imation, one has T>(t)/T>(0) = F(tN 2 ^ 3 ), with F(x) de- 
caying faster than l/x for large x, making the integral 
in equation converge. One has accordingly in the low 
T phase Jn — foa oc 1/iV 2 / 3 . A temperature indepen- 
dent exponent 2/3 for the free energy is in contradiction 
with the claims of ,13] that the internal energy scales like 
ejv — Coo oc 1/N xl - T \ with an exponent x(T) that is com- 
patible with 2/3 for both T = and T c but reaches a 
minimum w 0.54 between. The results of reference filjj 
are based however on Monte Carlo simulations of rela- 
tively small systems with N up to 196. Analyzing the 
data for the internal energy produced during the simu- 
lation of reference jl5jj_witch include systems with up to 
4096 spins, one finds 16] an exponent that is much closer 
to 2/3, with deviations that are presumably explained by 
the proximity of the critical point and by the very slow 
convergence of the expansion of e n — Soo in inverse powers 
of 1/N (at T c , the expansion parameter is[| I/A 1 / 3 ). 

The situation is different at T c , as shown in figure |3 
the ratio T>(t)/D(t — 0) scales with a different exponent, 
namely like F^N 1 / 3 ), with a large x behavior compat- 
ible with F(x) oc l/x (Although much larger system 
sizes would be needed in order to be sure that the sys- 
tem really approaches this asymptotic behavior). This 
is in agreement with formula 0] (In this model one has 
(3/v = 2). The data presented at T c (Figures |3J and [5J in- 
clude the results of an additional simulation of a system 
with N = 2048 sites, limited to the (cheap to simulate) 
paramagnetic phase, with T e [1.0,1.3], AT = 0.025, 
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FIG. 2: (Color online) V{t)/D{t = 0) as a function of £7V 2/3 
(both in logarithmic scale), for T = 0.6. The orange line (full 
line) shows the 1/x behavior, the blue line (dotted line) shows 
the 1/x 3 / 2 behavior. Clearly D(t) grows faster than 1/x for 
large x. The precise behavior of T>(t) is not essential for my 
argument, as soon as it decays faster than 1/x. 



with 128 disorder samples, and a 15 points discretization 
of t. 
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FIG. 3: (Color online) V{t)/D{t = 0) as a function of tN 1/3 
(both in logarithmic scale), for T = T c . The orange line 
(straight line) shows the expected 1/x behavior, in order to 
guide the eyes. 

Figure 0] shows, as a function of 1/iV 2 / 3 , my estimates, 
after integrating numerically equation^ of {fN~fN/2)/T 
at T = 0.4, compared to the result of a linear fit (/jv — 
fN/2)/T = -A/N- 2 ^, with A = 0.82 ± 0.02 and x 2 = 
4.9. The agreement is good within estimated statistical 
errors. A similar agreement is obtained for other values 
of T in the spin glass phase (e.g. A — 0.39 ± 0.01 with 
X 2 = 3.6 for T= 0.6, and A = 0.18 ±0.01 with X 2 = 33 - 
a large value presumably related to the proximity of the 
critical point - for T — 0.8). Figure [S] shows my estimates 










-0.005 




-0.01 






Z 


-0.015 


z 


-0.02 








-0.025 




-0.03 




-0.035 




0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 

1/N 2/3 

FIG. 4: (Color online) Numerical data for (/jv —/n/2)/T as a 
function of 1/N 2 / 3 , together with a numerical fit to the data 
of the form (f N - f N/2 )/T = -A/N 2/: \ with A = 0.82 ±0.02 
(blue dotted line). Here T = 0.4, N = 128, 256, 512 and 1024. 



f° r {In — In/2)/T at T c as a function of 1/N, together 
with the prediction of equation 01 A good agreement 
(with x 2 — 4.3 if one excludes the N = 128 data from 
the fit) is obtained using the value 1/N = 7.8 ± 0.2, 
namely /(_i) = ln(7.8)/12 = 0.17 . . ., within the range of 
results presented by Parisi, Ritort and Slanina 0. 
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FIG. 5: (Color online) Numerical data for (/jv — }n/2)/T as 
a function of 1/iV, together with the behavior implied by the 
equation: f N /T = f oa /T+ l/(12jV) laN/No. The orange line 
(full line) is drawn with the value No = 1. The blue line 
(dotted line) is drawn with the value 1/7.8, from a fit to the 
data. Here T = 1, TV = 128, 256, . . . , 2048. 

In conclusion, I have shown that the Guerra-Toninelli 
interpolation provides an efficient method to evaluate nu- 
merically the finite size corrections to the free energy of 
the Sherrington-Kirkpatrick model. The integrand T>(t) 
exhibits a remarkable scaling as a function of the inter- 
polation parameter t and system size N. At the crit- 
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ical temperature, the results for the free energy are in 
agreement with the predicted (1/12AT) \n(N/N ) leading 
behavior of the finite size corrections, and give the es- 
timate N w 1/7.8. In the low temperature phase, the 
results indicate that the leading corrections behave like 
TV -2 / 3 for both the internal energy and the free energy 
of the model. 
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